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ABSTRACT 

In  the  main  body  of  our  report, we  developed  an  asymptotic 
ray  theory  that  predicts  correctly  the  exterior  resonant  frequencies 
of  a  sphere.  A  logical  extension  of  this  work  is  to  spheroids 
where,  unlike  the  sphere,  the  polar  ray  orbit  resonance  condition 
is  nontrivial.  The  idea  is  thus  to  extend  the  results  to  separable 
coordinate  systems,  so  that  the  WKB  resonance  condition  can  be 
applied  independently  to  the  rays  along  the  three  coordinate  paths. 
The  radial  ray  path  remains  trapped  between  a  scattering  surface 
and  the  turning  point  caustic. 

At  this  point  the  results  of  this  addendum  are  only  formal 
and  tentative.  No  numerical  work  has  yet  been  carried  out. 


iv 


INTRODUCTION 


As  a  further  development  of  our  Ideas,  let  us  set  up  the 
exterior  resonance  conditions  on  a  prolate  spheroid.  Since 
natural  frequencies  are  available  in  the  literature  for  this  body, 
this  computation  can  be  checked  against  existing  theory  [1]. 

Thus  we  introduce  the  spheroidal  coordinates  (u,  9 , <f>  )  by 
stating  their  relation  to  the  Cartesian  coordinates  (x,  y,  z). 


x  ■  a  sinhu  sin  6  cos  <J> 

y  ■  a  sinhu  sin  6  sin  <p 

z  ■  a  coshu  cos  0 


(A.l) 


The  u  *  constant  surfaces  are  confocal  prolate  spheroids  rota- 
tionally  symmetric  about  the  z  axis.  In  the  limit  u  -*■  ®,  and 
a  •+■  0  such  that  r  ■  a  sinhu  is  finite,  the  system  (A.l)  returns 
to  the  usual  spherical  system  (r,  0,  <j>) .  Let  us  define  the 
contravarlant  vectors 


coshu 
cos  8 

4> 


(A.  2) 


Then  the  fundamental  arc  length  formula  is  defined  to  be 


ds 


Z  g..  du1,  du^ 


ij-l 


ij 


(A.  3) 


1 


where 


g 


ij 


.  I 

L  i  i 
4-1  3a  3uJ 


■  h  i  5iJ 

«2».V  -  (.Vl 

(u1)2  -  1 


2  a2 ((u2)2  -  (u2)2) 

2  1  -  <uV 


h23  -  a2((uV  -  1)  (1  -  (u2)2) 


Since  the  metric  tensor  g^  Is  diagonal,  the  prolate  spheroidal 
coordinate  system  is  seen  to  be  orthogonal.  From  the  scale  factors 
as  given  by  (A.  3)  one  computes  the  Laplacian  operator  according 
to  the  prescription 


V2  ij; 


JL  s  -i. 

/T  l-l  3u* 


3u 


(A. 5) 


where 


qi  -  TT*  .  g  -  h2  h2  h2  -  det(g±j) 
“  i 

2  2 

The  scalar  wave  equation  (V  +  k  )ip  -  0 

takes  the  form 


.  Li  1  -  2,2  '  A  K*"1)2  -  1,1  +  72  ld-(u2)2)  fL]} 

(u  )  “  (u  )  3u  3u  3u 


((u1)2  -  l)(l-(u3)2  3(u3)2 


32i)  2  (A.  6) 

(ka)2  ip  -  0, 
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Using  the  usual  separation  of  variables  techniques  we  let 


V  -  R(u1)  0  (u2)  $  (u3).  (A. 7) 

Substitution  of  (A. 7)  into  (A. 6)  eventually  leads  to  the  following 
three  separated  ordinary  differential  equations  for  R,  0  and  $. 

2 

«(uV  -  1)  rV1))'  +  ((ka)2  (u1)2  -  — S-a - +  A)R(uV  0 

(u’r-i 

(1  -  (U2)2)  0  '(u2))'  -  ((ka)2  (u2)2  + - X)  0  (u2)  -  0 

l-(u  )2 

®"  (u3)  +  m2  <f  (u3)  =  0  (A.  8) 

2 

In  (A. 8),  m  and  X  are  the  separation  constants.  The  equation 
3 

for  <f(u  )  requires  that  m  be  an  integer  to  satisfy  physical  constraint 

3  3 

that  u  *  0  be  the  same  point  as  u  3  2n. 

The  rotational  symmetry  of  the  prolate  spheroid  about  the  z 
axis  as  follows  from  equation  (A.l)  determines  the  separation  con¬ 
stant  m.  Consequently,  the  two  remaining  unknowns  for  the  geometric 
computation  of  resonant  frequencies  are  the  wave  number  k  and  the 
separation  constant  X. 

The  two  conditions  needed  to  determine  k  and  X  are  the  polar 
resonance  orbit  and  the  radial  resonance  orbit.  In  the  latter  con¬ 
dition  the  "orbit"  is  between  the  turning  point  caustic  and  the 
surface  of  the  spheroid. 


POLAR  ORBIT  RESONANCE  CONDITION 


Because  of  the  analytic  nature  of  the  phase  function,  there 
is  no  loss  of  generality  to  simplify  the  polar  resonance  condition 
by  computing  it  for  m  »  0  and  also  just  inside  the  caustic  surface. 

By  inspection  of  the  radial  differential  equation  for  R  given  in 
(A. 8),  the  turning  point  caustic  satisfies  the  equation 

X  ■  -  ft2  cosh2(u)  ,  ft  ■  ka  (A. 9) 

where  from  (A.l)  u  is  seen  to  be  the  radial  coordinate. 

In  (A. 9)  we  have  introduced  ft,  the  large  parameter  of  this 
asymptotic  calculation.  The  differential  equation  for  0  as  given  by 
(A. 8)  under  these  conditions  becomes 

((1-x2)  T' (x) ) '  +  ft2  (cosh2  u  -  x2)  T(x)  -  0  (A.10) 

9 

where  x  ■  u  =  cos  9  and  T(x)  »  0  (cos  9).  To  utilize  the  WKB 
approach,  we  make  the  further  change  of  variables 

y  -  tanh_1(x)  ,  ip  (y)  -  T  (x)  (A.  11) 

in  (A. 10)  to  obtain 

r  (y)  +  ft2  f2(y)  <y)  -  0  (A.  12) 

where 


4 


(A.  13) 


2  2 

f(y)  .  ^cosh  u  -  tanh  y^l/2 
cosh2y 

From  (A. 13)  it  is  seen  that  there  are  no  finite  turning  points.  Thus 
the  WKB  solution  of  (A. 12)  is 


$  (y)  *  exp  (i  QS(y)) 


(A.  14) 


where  the  phase  function  S(y)  is 

S(y)  =* 


f(y')dy’  +  ^  In  (f(y))  +  0(-y) 
y 


(A. 15) 


Returning  to  the  original  angle  variable  9  through  the  trans¬ 
formation  cos  9  =  tanhy,  the  real  (SR(0)),  and  imaginary  (5^.(9)),  parts 
of  the  phase  function  S(0)  become 

0 


SR(0) 


(cosh2u  -  cos29')^2  d9' 

0 


Si(9)  *  In  [(cosh2u  -  cos^)^2  sin  9  ] 


(A. 16) 


(A. 17) 


Let  us  interpret  this  result  geometrically  by  computing  a  geometric 
optics  solution  of  a  ray  traveling  along  a  polar  orbit  on  the  spheroid. 
The  geometric  optics  ray  solution  _  (0)  takes  the  form 

i  kl(e) 


*G.O.  (9)  *  e" 


where  1(9)  ■ 


9 


ds 


ds 


/ 


(dx)2  +  (dz)2 


(A.  18) 

(A.  19) 


From  (A.l)  the  arc  length  increment  ds  is  easily  seen  to  be 


2  2  1/2 
ds  ■  a(cosh  u  -  cos^  9)  d9 


(A. 20) 
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Thus  from  (A. 16)  and  (A. 18)  through  (A. 20) 


^G>o<(0)  =*  exp  (i  ft  SR  (0)) 


(A.  21) 


Comparing  the  geometric  optics  field  (A. 21)  with  the  WKB  solution 
(A. 14),  it  is  seen  they  differ  by  terms  of  order  0  (1/ft).  The 
attenuation  factor  S^(9)  is  of  order  0  (1/ft)  and  accounts  for  ray  tube 
spreading  on  the  curved  convex  ray  path. 

The  polar  resonance  condition  consists  of  requiring  that  the 
polar  phase  S(9)  satisfy  the  global  condition 


ft 


r  2ir 

S(9)  d8 

0 


2V 


n 


9 


1,2,3 - 


Because  of  four-fold  symmetry  this  can  also  be  written  as 

f-rr/2 


ft 


Jo 


S(9)  dG  =>  n.  tt/2 

U 


(A.  22) 


It  remains  to  evaluate  the  integral  in  (A. 22).  Separating  the 
phase  function  into  its  real  and  imaginary  parts  according  to  (A. 16) 
and  (A. 17)  produces  two  integrals.  The  imaginary  part  can  be  evaluated 
in  terms  of  elementary  functions.  One  determines 

tt/2 

SI(9)  d6  -  ~  [u  -  2  In  2]  (A. 23) 


The  corresponding  integral  of  S  (0)  is  expressible  in  terms  of 

R 

elliptic  integrals  of  the  second  kind.  This  is  not  surprising.  Elliptic 


6 


integrals  were  originally  encountered  in  the  calculation  of  the 
perimeter  of  an  ellipse  [4]!  Manipulation  of  (A. 16)  results  in 


S  (9)  »  csca  (E(sin  a) 
R 


-  E  (tt/2  -  6\a)) 


(A.  24) 


where 


a  -  Sin"1  ( — -r— ) 
coshu 


and  the  elliptic  integrals  of  the  second  kind  are  defined  as  [5]: 


E(<t>\a) 


2  2  1/2 
(1  -  sin  a  sin  9)  '  d9 


(A.  25) 


and  the  complete  integral  E(sin  a)  is  defined  as 

fff/2 


E(sin\a)  *  E(tt/2  a) 


(1  -  sin2  a  sin:  0)1/2d0  (A. 26) 


Substituting  (A. 24)  into  the  resonance  condition  (A. 22)  requires  the 
integral 


r*/2  „  Ft 2 

S  (0)d9  *  csca  {tt/2  ECsin^a)  -  0(1  -  sin^acosz0)  d9} 

o  Jo 

*  csca  {  [  (tt/2  -  9)  (1  -  sin2  a  cos2©)1^2  d0) 

'0 


(A. 27) 


Combining  results  (A. 22),  (A. 23)  and  (A. 27)  yields  the  resonance 
condition 


f^/2  2  2  1/2  iTT 

ft  csca  (  j  (tt/2  —  0 )  (1  —  sin^  acos  9  )  d0)  +  -£-(u-2  £n 2)  » 


where 


31x1-1  ('^kr)  » C08h2u 


\ 

7 


(A. 28) 
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Equation  (A. 28)  contains  two  unknowns,  X  and  Q.  The  required  addi¬ 
tional  condition  is  now  considered • 

RADIAL  RESONANCE  CONDITION 

As  was  described  in  the  main  body  of  the  report,  the  radial 
resonance  condition  requires  that  the  evanescent  ray  path  between  the 
scattering  surface  and  the  turning  point  caustic  be  an  integer  number 
of  wavelengths.  Included  in  this  phase  resonance  are  the  discontinuous 
phase  change  at  the  caustic  surface  and,  if  necessary,  at  the  surface 
of  the  scatterer. 

To  begin,  we  consider  the  radial  differential  equation  for 
RCu1)  as  given  in  (A. 8).  In  it  let 

s  =  coth  *  (u1)  ,  ip  (s)  •  R  (u1)  (A.  29) 

This  is  the  legimitate  change  since  u1  =  coshu  >  1.  The  equation  for 
^(s)  becomes 


lT(s)  +  s2  n2(ft)  <J>(s)  -  0 


(A.  30) 


where 

n2(s) 


coth2(s)  m2  |  \ 
sinh2(s)  n2 


2 

csch  s 


(A. 31) 


The  caustic  surface  is  defined  by  the  positive  value  of  s  denoted  by 

s„  such  that 
c 


n  (s  )  ■  0 
c 


(A. 32) 


The  solution  to  (A. 32)  is  the  proper  root  of 


8 


s  .  9illh-l  tSLti.--  if2*  V2  *  4fiV)1/2,l/2 

C  2a2 


m  >  0 


.  .-1  r  -n  , 
s  -  sinh  [— r - ] 

JT  +  X 


m  ■  0 


When  m  ■  0,  it  is  therefore  necessary  that  -X/Jr  >  1. 

For  eigenvalues  which  correspond  to  bound  discrete  modes  to 

exist,  it  is  required  that  the  effective  refractive  index  profile 
2 

n  (s)  have  the  property  that 


n  (s)  <  0 


s  <  s  <  s, 
c  b 


n  (s)  =■  0 


s  =  s 


(A.  34) 


n  (s)  >  0 


s  <  s 


In  (A. 34)  s^  defines  the  spheroid  surface.  The  radial  res¬ 
onance  condition  for  a  perfectly  conducting  prolate  spheroid  with  the 
electric  field  normal  to  the  surface  is 


I  n(s)ds  -  2tt  (nr  -  1/4) 


(A. 35) 


nr  =  1,2,3... 


In  (A. 35)  s^  is  the  surface  of  the  spheroid  defined  by 
coth(Sk)  ■  cosh(u^)  where  u^  defines  the  spheroid  through  equation  (A.l) 
when  u  ■  u^.  An  explicit  form  of  the  radial  resonance  is 


fc  c2- 

JPv 


P(s2  +  X)  -  P2S)2)1/2  — — 


P/P  ■+•  1 


27ri(nr  -  1/4) 
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(A.  36) 


where 


P.  *  cosh  u  -  1  and 
b  p 

2 

p.  -  1/2  [((i  +  x/ft2)2  +  >1/2  -  a  +  x/n2) ] 

n 

The  simultaneous  solution  of  (A. 28)  and  (A. 35)  for  X  and  ft 
constitutes  the  formal  asymptotic  solution  to  the  exterior  resonant 
frequencies  of  a  perfectly  conducting  prolate  spheroid. 
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